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Abstract 

A qualitative analysis is presented for spatially flat, isotropic and homogeneous 
cosmologies derived from the string effective action when the combined effects of 
a dilaton, modulus, two-form potential and central charge deficit are included. 
The latter has significant effects on the qualitative dynamics. The analysis is 
also directly applicable to the anisotropic Bianchi type I cosmology. 
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1 Introduction 



There has been considerable interest recently in the cosmological implications of string 
theory. String theory introduces significant modifications to the standard, hot big 
bang model based on conventional Einstein gravity and early universe cosmology 
provides one of the few environments where the predictions of the theory can be 
quantitatively investigated. A study of string cosmologies is therefore well motivated. 

The evolution of the very early universe much below the string scale and for 
string coupling much smaller than unity, g s <C 1, is determined by ten-dimensional 
supergravity theories |], ||, [|]. All theories of this type contain a dilaton, a graviton 
and a two-form potential in the Neveu-Schwarz/Neveu-Schwarz (NS-NS) bosonic 
sector. If one considers a Kaluza-Klein compactification from ten dimensions onto 
an isotropic six-torus of radius e 13 , the effective action is given by 



S = J d^x^-ge-* [i? + (V$) 2 - 6 (V/3) 2 - ^H^H^ - 2A 



where R is the Ricci curvature of the spacetime with metric g^ u , g = detg^, the 
dilaton field, $, parametrizes the string coupling, g 2 s = e*, and H^ u \ = d^B^] 
is the field strength of the two-form potential, B^ v . The volume of the internal 
dimensions is parametrized by the modulus field, f3. The moduli fields arising from 
the compactification of the two-form on the internal dimensions have been neglected 
Q. The constant, A, is determined by the central charge deficit of the string theory. 
In principle, it may take arbitrary values if the string is coupled to an appropriate 
conformal field theory. Such a term may also have an origin in terms of the reduction 
of higher degree form-fields || . 

In four dimensions, the three-form field strength is dual to a one-form: 

H ^x = eV^V^a, (1.2) 

where e^ uXK is the covariantly constant four-form. In this dual formulation, the field 
equations can be derived from the action 



S = J d A xy/—ge 



R + (V$) 2 - 6 (V/3) 2 - ^e 2 * (Va) 2 - 2A 



1.3) 



where a is interpreted as a pseudo-scalar 'axion' field ||. It can be shown that the 
action ( |1.3| ) is invariant under a global SL(2, R) transformation on the dilaton and 
axion fields when A vanishes 0. The general Friedmann-Robert son- Walker (FRW) 
cosmologies derived from Eq. (O) with A = have been found by employing this 



symmetry J/J. However, the symmetry is broken when a stringy cosmological constant 
is present M and the general FRW solution is not known in this case. 

The purpose of the present paper is to determine the general structure of the 
phase space of solutions for the class of spatially flat, FRW string cosmologies derived 



from the effective action ( |1.3| ) when a cosmological constant is present. This is well 
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motivated from a theoretical point of view and is also relevant in light of recent high 
redshift observations that indicate a vacuum energy density may be dominating the 
large-scale dynamics of the universe at the present epoch |J. 

The paper is organized as follows. In Section 2, the field equations are presented 
as an autonomous system of ordinary differential equations (ODEs). The combined 
effects of the axion, modulus and dilaton fields are determined for a negative and 
positive central charge deficit in Sections 3 and 4, respectively. This extends previous 



qualitative analyses where one or more of these terms was neglected ||T0|, |TT|, [12L O. H 



A full stability analysis is performed for all cases by rewriting the field equations in 
terms of a set of compactified variables. We conclude in Section 5 with a discussion 
of the phase portraits. 



2 Cosmological Field Equations 

The spatially flat, FRW cosmological field equations derived from action ( |1.3| ) are 
given by 



2a - 2a(p - a 2 e 2lp+&a 


= 


(2.1) 


2(p-^ 2 - 3d 2 - 6$ 2 + hy 2 e 2 ^ a + 2A 


= 


(2.2) 




= 


(2.3) 


a + & (ip + 6a) 


= o, 


(2.4) 



where 

ip = ^-3a (2.5) 

defines the 'shifted' dilaton field, a = e a is the scale factor of the universe and a dot 
denotes differentiation with respect to cosmic time, t. The generalized Friedmann 
constraint equation is 

3d 2 - 2 + 6/3 2 + h 2 e Wa + 2A = 0. (2.6) 

A number of exact solutions to Eqs. ( j2.1|) -( PTB|) are known when one or more of 
the degrees of freedom are trivial. We now discuss those that represent the invariant 
sets of the full phase space of solutions. The 'dilaton-vacuum' solutions, where only 
the dilaton field is dynamically important, are given by 

a = a*\t\ 1/p± 
e * = e "'\t\ p± -\ (2.7) 

where p± = ±Vo and {a*, $*} are arbitrary constants. There is a curvature singular- 
ity at t = 0. The solution ( |2.7| ) forms the basis of the pre-big bang scenario, where 
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a growing string coupling can drive an epoch of inflationary (accelerated) expansion 
15 1 . The pre-big bang phase corresponds to the p = p_ solution over the range t < 
and the post-big bang phase to the p = p + solution for t > |T5]. The inflationary 
nature of this scenario has recently been questioned, however [0| . 

The 'dilaton-moduli-vacuum' solutions, with a = A = 0, are given by 



a* \t 



±h, 



$ r.|±3/i,-l 

e = e \t\ 



J = e P* ui±V(i-3ft2)/6 



(2.8) 



where {/t*, are constants. This class of solution can also be expressed in the form 

, (l+r±)/2 



s ^ 



r± 



(2.9) 



where {s*, g, r±} are constants and s = /* dt'/a(t') is conformal time. The generalized 
Friedmann constraint equation (|2.6|) leads to the constraint r± = ± (3 — ^g 2 ) 1 ^ 2 . 
The general solution where only A = is the 'dilaton-moduli-axion' solution [[7]: 

1/2 
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1/2 
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—r 
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e * 
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(x* ± e 



+ 
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5^ 










— r 


s/s* \ r 




Js/ S * 


-r + 


|s/s*| r 



(2.10) 



where <t* is an arbitrary constant and r = |r±|. This cosmology asymptotically 
approaches one of the dilaton-moduli-vacuum solutions ( |2.9| ) in the limits of high 
and low spacetime curvature. The axion field induces a smooth transition between 
these two power-law solutions and causes a bounce to occur. It is only dynamically 
important for a short time interval when s ~ 0. 

The solutions where only the axion field is trivial and A > are specific cases of 
the 'rolling radii' solutions [ |I7|j : 



a 



a* |tanh(At/2)| r 



-<!>,, 



cosh(At/2)| 2fc - 6n |sinh(At/2)| 



2Z+6n 



|tanh(At/2)| ? 



(2.11) 
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where A = v2A and the real numbers {k, I, m, n} satisfy the constraints 



3m 2 + 6n 2 = 1, 3m + 6n = k — L k + I = 1. 



(2.12) 



The corresponding solutions for A < are related to Eqs. ( 2.11| ) by redefining A 
—iA. In this case, the range of t is bounded such that < t < tt/A. 



Finally, there exists the 'linear dilaton vacuum' solution where A > |T8] . This 
solution is static and the dilaton evolves linearly with time: 



0. 



p = o, 



$ = ±V2At. 



(2.13) 



The field equations 
tonomous ODEs: 



T|)- (|2.6|) may be written as the following system of au- 



h 

N 
P 
3h 2 



= ip 2 + htp - 3h 2 - N - 2A 

= 3h 2 + N 

= 2Nip 

= —6hp 

- tp 2 + N+ lp + 2A = 0, 



where we have defined new variables 



N = 6P 2 



i[) = cp, 



h 



a. 



(2.14) 
(2.15) 
(2.16) 
(2.17) 

(2.18) 



(2.19) 



The variable p may be interpreted as the effective energy density of the pseudo- 
scalar axion field [ 13| . It follows from Eq. ( |2.15|) that ip is a monotonically increasing 
function of time and this implies that the equilibrium points of the system of ODEs 
must be located either at zero or infinite values of if). In addition, due to the existence 
of a monotone function, it follows that there are no periodic or recurrent orbits in the 

The sets A = and p = are invariant sets. In particular, the 
given by Eqs. ( |2.1U| ) divides the phase space and the orbits 



phase space 19, 20 



exact solution for A 

do not cross from positive to negative A. 

We now proceed to consider the cases where A < and A > separately. 



3 Analysis for Negative Central Charge Deficit 
3.1 Four— dimensional Model 

In this Section we consider the phase portraits of the NS-NS fields for negative central 
charge deficit. In the case where the modulus field is frozen, = 0, Eq. (|2.18|) may 
be employed to eliminate the axion field's energy density. This reduces the set of Eqs. 
fl2.14j )- (|2.171 ) to a two-dimensional system. Moreover, it follows from Eq. ( |2.18|) that 

^ 2 -2A>3/i 2 >0 (3.1) 
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and we may therefore compactify the phase space by defining new variables 



V = 



yj$* - 2A 

V^ 2 - 2A 



and a new time variable 



d _ 1 d 
Eqs. ( 2.14 ) and ( 2.15 ) then become: 



g??7 

d£ 
dr 



1 - 



(3.2) 
(3.3) 

(3.4) 

(3.5) 
(3.6) 



The variables defined in ( 3.2 
it follows from Eqs 



and (|3.3|) are bounded, rj < 1 and £ < 1, and 
(|3.5| ) and (|3.6j ) that they are both monotonically increasing 
functions. The equilibrium points are located at £ 2 = rj 2 = 1. The invariant sets 
p = and A = correspond to the conditions £ 2 = 1 and rf — 1, respectively. A 
stability analysis indicates that the equilibrium point A : (77, £) = (1, 1) is an attractor 
and the point J? : (77, £) = (—1, —1) is a repeller. The points Si i2 : (1, —1) and (—1, 1) 
are both saddles. The phase portrait is given in Fig. 1 and is discussed in Section 5, 
where a physical interpretation is given. 



[FIGURE 1 HERE] 



3.2 Ten— dimensional Model 

We now consider the effect of lifting the solutions to ten dimensions by including the 
modulus field, (3. This will raise the dimension of the phase space to three. In this 
case, it proves convenient to employ the generalized Friedmann constraint equation 
( p. 181 ) to eliminate the modulus field rather than the axion field. This equation can 
be written as 

— 5^5a- (3 ' 7) 

where the new variable k is defined by 

« = ^-^r- — rr (3-8) 
2(ij 2 - 2A) v ; 
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and satisfies < k < 1. Employing Eqs. 
equations (|2.14|) - (|2.17|) in the form: 



dt] 

dr 

dr 
dn 



(1 - k) [1 - rf 



— = -2k 



v^ + rKl 



we can now express the field 

(3.9) 
(3.10) 

/.') . (3.11) 



We note that the invariant set N = corresponds to k = 1 — £ 2 , in which case the 
above system of ODEs reduces to the two-dimensional system (|3.5|) - (|3.6|) . 

The equilibrium points of this system of ODEs all lie on one of the two lines of 
non-isolated equilibrium points (or one-dimensional equilibrium sets) 



L± : 7)1 = 1,k = 0, 



(3.12) 



where £ is arbitrary. The corresponding eigenvalues are Ai = — 2r/ and A2 = 
— 2(\/3£ 4- rjo), and hence these equilibrium sets are normally hyperbolic (through- 
out, we shall refrain from giving the corresponding eigenvectors explicitly). The 
third eigenvalue is zero since this is a set of equilibrium points. Thus, on the line 
L + : (r)o — 1, k — 0; £) the equilibrium points are saddles for £ G [ — 1, — 1/ -y/3) and 
local sinks for £ G (— 1/Vo, 1]. On the line L_ : (r) = —1, k = 0; £) the equilibrium 
points are local sources for £ G [—1, l/\/3) and saddles for £ G (l/y/3, 1]. The phase 
portrait is given in Fig. 2. The dynamics is very simple due to the fact that the 
right-hand sides of Eqs. (|3.9| ) and (|3.10|) are positive-definite and hence 77 and £ are 
always monotonically increasing functions. The curved upper boundary k = 1 — £ 2 
denotes the invariant set N = and therefore corresponds to Fig. 1. 



[FIGURE 2 HERE: 'LARGE'] 



4 Analysis for Positive Central Charge Deficit 

In the case where the central charge deficit is positive, the variable tp 2 — 2A is no 
longer positive-definite and therefore can not be employed to normalize the system. 
In view of this, we choose the normalization 

/ 1 \ 1/2 

e = [3h 2 + -P + N + 2A) . (4.1) 
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The generalized Friedmann constraint equation Q2.18] ) now takes the simple form 

(4.2) 



ti 2 

= 1 



and may be employed to eliminate ip. Since by definition e > 0, specifying one of the 
roots ip/e = ±1 corresponds to choosing the sign of ip. However, it follows from the 
definition in Eq. ( |2.19|) that changing the sign of ip is related to a time reversal of the 
dynamics. In what follows, we shall consider the case ip/e = +1; the case ip/e = —1 
is qualitatively similar. 

Introducing new variables 



V3h 



v 



and a new dynamical variable 



d_ 

dT 



= _P_ 

~ 2e 2 ' 

1 d 



transforms Eqs. (|2.14|) - (|2.17|) to the three-dimensional autonomous system: 



d/j 
dT 
dv 
dT 
dX 
dT 



v + 



-2u 



V3 



1 — /i 
1 



1 



A 



(4.3) 
(4.4) 

(4.5) 
(4.6) 
(4.7) 



The phase space variables are bounded by the conditions < {^ 2 ,u, A}<1 and also 
satisfy fi 2 + v + A < 1. The sets v = and A = are invariant sets corresponding 
to p = and /3 = 0, respectively. In addition, fi 2 + v + A = 1 is an invariant set 
corresponding to A = 0. We note that the right-hand side of Eq. (|4.7|) is positive- 
definite and this simplifies the dynamics considerably. 

The equilibrium points of the system ( |4.5| )-( |4"7F| ) consist of the isolated equilibrium 
point 

C-.fx = u = X = (4.8) 
and the line of non-isolated equilibrium points 



V 



0,A 



jj 2 (/i arbitrary). 



(4.9) 



The eigenvalues associated with C are Ai = l/y/3, A 2 = 2/y/3 and A 3 = 0. Although 
this isolated singular point C is non-hyperbolic, a simple analysis shows that it is a 
global source. The eigenvalues associated with V are: 



Ai 




-=], A, = = 



2 



(4.10) 
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and the third eigenvalue is zero since V is an equilibrium set. Therefore, on V the 
equilibrium points are saddles for p, G [— 1 , — 1/ v^3) and local sinks for /x G ( — 1 / 1] . 
The phase portrait is given in Fig. 3. 



[FIGURE 3 HERE: 'LARGE'] 



It is also instructive to consider the dynamics on the boundary corresponding to 
A = 0, since the case N = is of physical interest in its own right as a four-dimensional 
model. In this case the ODEs reduce to the two-dimensional system: 



dT */a V p 



dp 
dT 



— 2/xz/ 



1 + 7T 



0: 



A, 







The equilibrium points and their corresponding eigenvalues are 

1 

73' 
2 

73 
2 

71 



C : 



S : 



-1,1/ 



Ai 
0; Ax 



A : /x = 1, v = 0; Ai 



~=, A, = 2 1--= 



A 5 




(4.11) 
(4.12) 

(4.13) 
(4.14) 
(4.15) 



Point C is a non- hyperbolic equilibrium point; however, by changing to polar coordi- 
nates we find that C is a repeller with an invariant ray 6 = tan _1 (— y/3). The saddle 
S and the attractor A lie on the line V. The phase portrait is given in Fig. 4. 



[FIGURE 4 HERE] 



This concludes the analysis of the phase portraits for the FRW string cosmologies 
containing non-trivial NS-NS fields. 

5 Discussion 

The phase portraits have a number of interesting features. In Fig. 1 the modulus 
field is frozen and the universe contracts from a singular initial state. The orbits 
in the vicinity of the equilibrium point R are asymptotic to the p = p_ dilaton- 
vacuum solution Q2.7| ). The axion is negligible and the kinetic energy of the dilaton 
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dominates the energy-momentum tensor. As the collapse proceeds, however, the 
axion becomes dynamically more important and eventually induces a bounce. In the 
case of vanishing A, Eqs. fl2.10|) imply that the future attractor would correspond 
to the p = p+ dilaton-vacuum solution. However, the combined effect of the axion 
and central charge is to cause the universe to evolve towards the equilibrium point 
A, where at — > +00, in a finite time. This behaviour differs from that found when no 
axion field is present, because in this latter case there is no bounce 0, IT7 . 



This behaviour can be understood by viewing the axion field as a membrane [21 



Since this field is constant on the surfaces of homogeneity, the field strength of the 
two-form potential must be directly proportional to the volume form of the three- 
space. If the spatial topology of the universe is that of an isotropic three-torus, the 
axion field can be formally interpreted as a membrane wrapped around this torus [ pTf . 
As the universe collapses, the membrane resists being squashed into a singular point 
and this forces the universe to bounce into an expansionary phase. The cosmological 
constant then dominates the axion field as the latter's energy density decreases. 

The inclusion of a modulus field leads to a line of sources and sinks for the orbits 
(see Fig. 2). The axion field is dynamically negligible in the neighbourhood of the 
equilibrium points. Moreover, a bouncing cosmology is no longer inevitable and there 
exist solutions that expand to infinity in a finite time. The solutions are asymptotic to 
the dilaton-moduli- vacuum solutions (|2.9|) near the lines L±. The boundary points 
£ 2 = 1/3 on these lines correspond to the limiting cases where a 2 = j3 2 . These 
represent the isotropic, ten-dimensional cosmology (a = $) and its dual solution 
(a = —(3). In the latter solution, the ten-dimensional dilaton field, $ = $ + 6/3, is 
constant. 

In Figs. 3 and 4, the isolated equilibrium point C corresponds to the 'linear 
dilaton vacuum' solution ( |2.13|) fill , HH- When the modulus is frozen, all trajectories 
evolve away from C towards the point A and approach the superinflationary p = p_ 
dilaton-vacuum solution ( |2.7| ) defined over t < 0. Some of the orbits evolving away 
from C represent contracting cosmologies and the effect of the axion is to reverse 
the collapse in all these cases. For the rolling modulus solution (Fig. 3), the orbits 
tend to the dilaton-moduli-vacuum solutions as they approach the attractors (the 
sinks on V). As in the case of a negative central charge, the critical value fi 2 = 1/3 
corresponds to the case where a 2 = f3 2 . The other boundary of V is the point A 
representing the case where the kinetic energy of the modulus field vanishes. The 
qualitative behaviour of models with ip < is similar. 

The results presented in this paper can be directly applied to the class of spatially 
flat and homogeneous Bianchi type I models by reinterpreting the physical meaning of 
the modulus field, j3. The line element for the axisymmetric type I model is given by 
ds 2 = —dt 2 + h a t,(t)dx a dx b (a, b = 1, 2, 3), where the metric on the surfaces of homo- 
geneity is defined by h a t, = e 2ct diag[— 2/5, (3, [3] f22j. The field equations derived from 



action (1.1) for this background are formally identical to those presented in Section 2 



In this anisotropic model, however, the variables a and (3 now parametrize the 
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averaged scale factor and the shear parameter of the universe, respectively. It would 
be interesting to investigate whether the above results can be employed to study the 
question of isotropization in string cosmologies. 

Finally, we observe that all of the exact solutions corresponding to the equilibrium 
points of the governing autonomous systems of ODEs are self-similar since in each 
case the scale factor is a power-law function of cosmic time (see, for example, Eq. 
PH. This result can be proven in general. For example, all of the equilibrium 



points of the system of ODEs ( |3 . 9| ) - (|3TTT| ) necessarily have h = constant (leading to a 



power-law solution) or A = 0, whence from the exact 'dilaton-moduli-axion' solution 
( p.lO| ) we can see that asymptotically the scale factor is power-law. This implies 
that exact self-similar solutions play an important role in determining the asymptotic 
behaviour of string cosmologies (19| . 

In conclusion, therefore, we have presented a qualitative analysis of the spatially 
flat, FRW string cosmologies containing non-trivial dilaton, axion and modulus fields 
together with a central charge deficit (stringy cosmological constant). In all cases, 
variables were found that led to a compactification of the phase space and this allowed 
a complete stability analysis to be performed. Both four- and ten-dimensional models 
were studied by including the dynamics of the modulus field. The combined effects of 
the axion field and central charge deficit on the qualitative behaviour of the dilaton- 
moduli-vacuum solutions (|2.8|) were determined. We found that such terms have a 
significant effect on the dynamics. 
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Figures 



Figure 1: The phase portrait of the system (|3.5|)-(|3^) corresponding to the four- 
dimensional NS-NS model with no modulus field and negative central charge deficit 
(A < 0). Equilibrium points are denoted by dots and the labels in all figures cor- 
respond to those equilibrium points (and hence the exact solutions they represent) 
discussed in the text. We shall adopt the convention throughout that large black 
dots represent sources (i.e., repellers), large grey- filled dots represent sinks (i.e., at- 
tractors), and small black dots represent saddles. Arrows on the trajectories have 
been suppressed since the direction of increasing time is clear using this notation. 



Figure 2: The phase portrait of the system ( p.9|) ( p7TT| ) . This corresponds to the 



ten-dimensional NS-NS model with negative central charge deficit (A < 0). Grey 
lines represent typical trajectories found within the two-dimensional invariant sets, 
dashed black lines are those trajectories along the intersection of the invariant sets, 
and solid black lines are typical trajectories within the full three-dimensional phase 
space. Note that L± denote lines of non-isolated equilibrium points. See also caption 
to Fig. 1. 



Figure 3: The phase portrait of the system (|4.5| )- (|4.7| ) corresponding to the ten- 
dimensional NS-NS model with positive central charge deficit (A > 0). The root 
ip/e = +1 of Eq. ( f4.2|) has been chosen. Note that V denotes a line of non-isolated 
equilibrium points. See captions to Figs. 1 and 2. 



Figure 4: Phase portrait of the system ( |4.11| )-( fl.l2j ) corresponding to the four- 
dimensional NS-NS model with positive central charge deficit (A > 0). The root 
tp/e = +1 of Eq. (|4.2|) has been chosen. See caption to Fig. 1. 
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NOTE TO EDITOR: Figs. 2 and 3 ARE TO BE AS "LARGE" AS POSSIBLE. Figs. 
1 and 4 can be "smaller" . 
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